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Two questions connected to the macroscopic Maxwell equa- 
tions are addressed: First, which form do they assume in the 
hydrodynamic regime, for low frequencies, strong dissipation 
and arbitrary field strengths. Second, what does this tell us 
about irreversibility and coarse-grained description. 
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I. INTRODUCTION 



You are probably ambivalent toward the article you 
have just started to read, or even harbor dark suspicions 
about it. These feelings are admittedly hard to avoid 
when encountering someone purporting to reveal news 
about the centenarian Maxwell equations. Afterall, we 
know all there is to know about the Maxwell equations, 
do we not? Well, no, we do not. Despite the impres- 
sion we came away with from various courses on elec- 
trodynamics, there are large gaps in our understanding 
on this subject, and the most glaring ones deal with the 
connection between electrodynamics and thermodynam- 
ics. These two independently developed classical areas 
of physics apparently do not mix well: Jackson (1975) 
does not mention entropy in his classic book at all, while 
Callen (1985) has eliminated the (partly erroneous) chap- 
ter on magnetic and electric system in the second edition 
of his definitive work. 

Nevertheless, it must be sensible to ask questions such 
as when an electromagnetic field configuration is in equi- 
librium, and how this equilibrium state is arrived at. 
Equilibrium fields are the ones that maximize the en- 
tropy, is there a simple equation expressing this property? 



Equilibrium fields are always static, but is it possible, in 
dielectrics, for a static field to be off-equilibrium? These 
are some of the questions answered below, and we go on 
from these answers to understand how field dissipation 
can be accounted for in dynamic situations, for a general 
system of nonlinear constitutive relations - where pos- 
tulating imaginary parts for the permeabilities e and fi 
fails to work. We also go on to understand what the elec- 
tromagnetic force is that a polarizable and magnetizable 
body feels - both in and off equilibrium. Fortunately, 
all these results are fairly simple and universal, and in 
fact quite suitable for being introduced into a university 
course on "electrodynamics in continua" . 

In section O, I shall summarize the state of the art 
of the Maxwell equations as it is conventionally treated, 
and point out the extent of securely understood physics. 
Starting with section HI , a different approach, that of the 
thermo- and hydrodynamic theory, is introduced. The 
hydrodynamic theory, the response of dense and dissi- 
pative systems exposed to slowly varying external fields 
of nonlinear strength, is discussed in section [v|. Inter- 
spersed between these considerations, I shall often stop to 
deliberate over two fundamental aspects of macroscopic 
theories — coarse-grained description and irreversibility 
— starting right at the beginning of the next section, but 
mainly in section VI. The difference between the hydro- 
dynami c an d the linear response theory is dwelt on in 
section VII, with some surprising and instructive results. 



Finally, in section VIII, the old but confusing subject 



about the electromagnetic force in a coarse-grained de- 
scription is considered and clarified. 

[Arguments that are not usually reproduced in a semi- 
nar, only invoked when the appropriate question is raised, 
are given - as here - in large square brackets. They may 
be skipped on first reading.] 



II. LESS ACCURATE IS MORE DIFFICULT? 

The actual difference between the electrodynamics in 
vacuum and continua is one of accuracy, better: one of 
resolution. Denoting the grain size of descriptive grid as 
£, and the average distance between the charge carriers 
as d, the vacuum Maxwell equations as a high resolution 
theory are valid for d ^> £■ Since the input, the density 
of charge and current, is in principle arbitrarily accurate, 
we are able to calculate the fields e and b to any desired 
resolution. 

Their contribution to the energy, and the Lorentz force 
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are, respectively 

u cm = i( £o e 2 + 6 2 /Vo), t = g e (e + vxh). (1) 

In conjunction with the Newtonian equation of motion, 
these two expressions account for the feedback, for how 
the field affects the motion of the material. So we have 
at our disposal a closed theory — a classical one with 
a known quantum mechanical generalization — that is 
conceptually simple yet technically intractable for dense 
systems. 

The circumstances are reversed if we check our ambi- 
tion and seek enlightment from a low-resolution theory, 
d <C £: The technical difficulties are greatly reduced, but 
conceptually we enter murky waters. 

Because the density of charge (g e ) and current (j e ) = 
(g e v), are now spatially averaged quantities, we must deal 
with hidden charges and currents, making it necessary to 
consider four instead of two fields: D, B, E M , and H M 
- all coarse-grained, hence in capital letters. (As will 
become clear soon, we need to distinguish E M , H M from 
E, H. The former two are the usual fields as defined by 
the Maxwell equations; the latter two will be introduced 
later.) Given constitutive relations, linear if the field is 
sufficiently weak, 

D = ( £ ' + * £ ")E M , B = (// + i M ")H M , (2) 

we are again able to calculate the field from the source. 
The real parts of the permeabilities, s' and p! , account 
for the reactive responses such as the oscillatory mo- 
tion of hidden charges; the imaginary parts, e" and 
fx", parametrize dissipation and absorption. If a field is 
stronger, it is customary to take the corresponding per- 
meability again as a function of the field. But one goes 
beyond linear response only at the price of loosing all 
the simple relations, especially the identification of the 
imaginary part with dissipation. 

Deplorable as this is, the problems are worse for the 
feedback, the effect of the field on the motion of the ma- 
terial. The reason is that both the electromagnetic en- 
ergy and the Lorentz force, Eqs ([!]), are nonlinear. And 
the knowledge of the coarse-grained quantities (g e ), (j e ), 
E M = (e) , B = (b) ... is quite useless if we need to know 
the values of (g e e), or (g e v x b). Nevertheless, bold souls, 
without much ado, simply write 

p&m + W t P + VjTrg = [(g e )E M + (je) x B] . . (3) 

In the absence of fields, the right side of Eq(||) is zero, and 
what remains is the Navier-Stokes equation, an expres- 
sion of momentum conservation in the low-resolution, hy- 
drodynamic physics: Both the pressure P, a quintessen- 
tially thermodynamic quantity, and the viscous stress 
tensor it®, a dissipative term, presume an infinitesimal 
volume element that contain enough particles to form a 
system in local equilibrium. This volume element is noth- 
ing but the descriptive grain, of the size £ 3 , we therefore 
have d <C £. Now, any theory, and each equation, must 



consistently have a unique resolution. So the sources and 
fields on the right hand side must not be of high reso- 
lution, but it is not at all clear that the expression as 
written, in low-resolution quantities, is correct. Needless 
to say, anyone employing this fairly popular equation, or 
a variant of it, bears the burden of proof for its validity. 
[An example of a system in which Eq(||) does hold is a 
weakly dissociated gas, where the density of charge car- 
riers is much lower than the density of neutral particles. 
Assuming negligible dipole moments for the latter, we 
have two interparticle distances, d n for neutral particles, 
and d e for the charge carriers. Then the resolution of 
this equation may be chosen as d n <C £ -C d e , such that 
it is of low-resolution for the neutral particles, yet of high 
resolution for the charge carriers.] 

Actually, things are not as bleak as they seem right 
now, and a few lesser known results do considerably 
brighten up the prospect of formulating the feedback, 
and closing the low-resolution Maxwell theory. These re- 
sults are (i) the thermodynamics of the electromagnetic 
field and (ii) the expressions for the energy and stress 
tensor at finite frequencies, albeit without dissipation. 
They will be briefly outlined in the next section and are 
taken from the one good book on this subject, Volume 
VIII of Landau and Lifshitz (1984), referred to below as 
LL8; see also Kentwell and Jones (1987) for a nice review 
of the state of the art. The thermodynamic considera- 
tions as presented, however, contain enough generaliza- 
tion and shifts in interpretation, from the unfortunately 
brief treatment in §18 of LL8, that I most probably have 
to bear the blame if you should find faults in any state- 
ments here. 



III. THE THERMODYNAMICS OF FIELDS 

We start from the energy density u in the system's 
rest frame, as a function of the entropy density s, mass 
density p, and the two fields D and B, 

du = Tds + fidp + E • dD + H • dB. (4) 

Being defined as du/dT) and <9w/<9B, respectively, E and 
H are (like T = du/ds) functions of all the thermody- 
namic variables. For weak fields (and barring ferroelec- 
tricity or ferromagnetism) , an expansion yields 

E = D/e, H = B//2, (5) 

where e and p may be functions of temperature and den- 
sity, but not the frequency. 

Maximizing the entropy S — JdVs for a stationary di- 
electric medium, with the constraints of constant energy, 
constant mass, and the validity of the two non-temporal 
Maxwell equations, V • D = (g e ) and V • B = 0, the 
resulting Euler equations are 

VT = 0, Vfi= 0, (6) 
V x E = 0, V x H = 0. (7) 
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The entropy is maximal, and the system in equilibrium, 
only if these equations are satisfied. [As an illustration, 
keep only D as a variable and minimize U — J u. This 
leads to /E-JD - $(r)<5(V-D - (g e )) = 0, with $(r) 
a Lagrange parameter. A partial integration and the 
fact that <5D is arbitrary yield E = — V<fr, or the first 
of Eqs(0). Circumstances are modified if a small but fi- 
nite conductivity eventually allows the charge to move, 
rendering (g e ) no longer constant, only J (g e ). Then the 
entropy can be further increased, and becomes maximal 
for E = 0.] 

The two Eqs(0) are a surprise, as they amount to a 
thermodynamic derivation of the static Maxwell equa- 
tions. Irrespective of what specific form u — and hence 
E and H — assume, Eqs(Q) hold. [To see that E and 
H are indeed the measured fields, we remind ourselves 
that the two variables D and B must reduce to e and b 
in vacuum, or a rarefied gas. A comparison of the two 
energy expressions, Eq(Q) and (0) then compels E and H 
to do the same, reduce to e and b. The usual boundary 
conditions then imply that these four fields of a dense 
system are the ones that will be appropriately continued 
into an adjacent vacuum, where field measurements are 
easily carried out.] 

This thermodynamic consideration - in conjunction 
with its hydrodynamic generalization in section |V] - pries 
open a new door to understand the macroscopic Maxwell 
equations. It possesses all the advantages of a thermo- 
and hydrodynamic theory, being general, independent 
of microscopic interactions, and valid for arbitrary field 
strength. At the same time, it is irreversible, and intrin- 
sically of low resolution. 

In the presence of fields, particles can be moved and 
accelerated, and the momentum density pv of the mate- 
rial is not a conserved quantity. However, because the 
Hamiltonian remains invariant under translations that 
include both the material and the charges producing the 
field, the total momentum of material and field is still 
conserved, 



v,n,, = 0. 



This total momentum density is 

g tot = p« + Ex H/c 2 , 



(8) 



(9) 



because (i) the main contribution in the relativistic, total 
energy current is the sum of the rest mass motion and the 
Poynting vector, Q rcl = [pc 2 +0(u)]t>+ExH; and (ii) the 
4-energy-momentum tensor is symmetric, Q rcl /c 2 = g tot . 
Despite controversies about the form of g tot that refuses 
to die down, it appears difficult to circumvent this simple 
and fundamental argument. 

On the other hand, for the gist of this paper, it is not 
important which form g tot assumes, as long as it is a def- 
inite one. The point is, given the expressions for g\ ot and 
Hij, the acceleration -^(pvi) may be calculated. There- 
fore, the knowledge of these two expressions contains that 
about the coarse-grained Lorentz force. 



The flux Hij is, in equilibrium and at vanishing veloc- 
ity, the Maxwell tensor. It may be derived by consider- 
ing the change in the total energy, of a system containing 
charges, when a certain portion of its surface is moved, 
see §15 of LL8. The result is a longish expression con- 
taining only thermodynamic and conjugate variables, 



U i:j = (Ts + [ip + E i D l + HiBi - u)5ij 
- [E.Dj + HiBj j)}/2. 



(10) 



The expression for the flux at finite frequencies is a 
less trivial matter and presupposes the corresponding 
expression for the energy. Assuming linear constitu- 
tive relations, lack of dissipation (e" , p" = 0), quasi- 
monochromacy and stationarity (v = 0), Brillouin 
showed in 1921 that the additional energy due to the 
presence of fields is 



±(E 2 )d(ujE')/dLu + i(# 2 )d(w//)/dw, 



(11) 



where the average is temporal, over a period of oscilla- 
tion. Compared to the corresponding thermodynamic ex- 
pressions, EqsQ |5|), there are two new terms ~ wde'/dw, 
wdp' /du. Forty years later, Pitaevskii showed that un- 
der essentially the same assumptions, the stress tensor 
retains its form from equilibrium, Eq(|Io|), and remark- 
ably, does not contain any frequency derivatives, cf §80, 
81 of LL8. 



IV. SOME NEW RESULTS 

If we draw a diagram of field strength versus frequency, 
with the field strength pointing to the right, and u> point- 



ing upward, see Fig 1, 
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FIG. 1. Range of validity 
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we have a vertical stripe A along the w-axis — arbi- 
trary frequency but small field strength — that is the 
range of validity of the linear response theory, Eq(||), 
while the field-axis itself depicts the space in which 
thermodynamics holds; the expressions of Brillouin and 
Pitaevskii are valid within the linear response stripe A, 
in isolated patches disjunct from the field-axis, wherever 
field dissipation is negligible. [Because s",fj," are odd 
functions of u, and e', \j! even ones, and because the first 
frequency dependent effects when leaving the equilibrium 
are linear in uj, they belong to e" , /i" and are dissipative. 
Lack of dissipation therefore characterizes a frequency 
region disjunct from uj = 0.] 

Having understood the thermodynamic behavior of a 
system, it is fairly easy to derive the corresponding hy- 
drodynamic theory. It accounts for the same system — 
condensed matter, charged or exposed to an external field 
— that is now slightly out of equilibrium, to linear order 
in the frequency. Its range of validity is the horizontal 
stripe B along the field axis, for arbitrary field strength 
and small frequencies. The theory as derived (Liu 1993, 
1994, 1995) is closed; it includes both the macroscopic 
Maxwell equations and the expression for the total mo- 
mentum flux Hij- As terms of second order in the fre- 
quency are neglected, the hydrodynamic theory does not 
account for dispersion. 

The parameter space C beyond the two perpendicular 
stripes needs a theory that can simultaneously account 
for dissipation, dispersion, nonlinear constitutive rela- 
tions and finite velocities. Although one might expect 
principal difficulties in setting up such a theory — due 
to the apparent lack of a small parameter — the system 
is in fact, up to the optical frequency ~ 10 15 Hz, still in 
the realm of macroscopic physics, as the associated wave- 
length remains large compared to the atomic graininess. 
And when asking questions such as what is the force on 
a volume element exerted by a strong laser beam, if we 
confine our curiosity to the averaged force — with a tem- 
poral resolution larger than the time needed to establish 
local equilibrium (which itself is much larger than the 
light's oscillatory period for the frequency range under 
consideration) — a simple, universal and hydrodynamic- 
type theory is still possible. Better: it may be cogently 
derived, as the respective limits of small field and low 
frequency are firmly anchored. A first step toward such 
a theory has been successful. It includes the dynamics of 
polarization, but neglects magnetization (Jiang and Liu 
1996). 

V. THE HYDRODYNAMICS OF FIELDS 

Although the hydrodynamic theory of electromag- 
netism may appear unorthodox at times, it is rather 
elementary in its essence, and especially easy to com- 
prehend by analogy. Consider a typical hydrodynamic 
equation, that of Navier-Stokes in the absence of fields, 



9i + ^ j^ij = 0. The momentum density gi = pvi is a 
thermodynamic variable, odd under time inversion. The 
stress tensor Ily = mj + irB is the corresponding flux, 
with two parts: The reactive one is (if linearized) given 
by the pressure, 7r^ = PSij, a thermodynamic derivative. 
It is even under time reversal, same as <?;. The dissipative 
part, 

is odd and breaks the time inversion symmetry of the 
equation. If we had included the momentum density gi 
as an additional variable in the thermodynamic consider- 
ation above, v^j = would have been added to Eqs(|[0) 
as the respective Euler equation. If it is satisfied, the 
entropy is maximal with respect to the distribution of gi. 
So understandably, if v^j ^ 0, there is a current ~ Uy 
to redistribute gi, such that the system is pushed toward 
maximal entropy and equilibrium. This is how dissipa- 
tion and irreversibility are generally accounted for in hy- 
drodynamic theories. Every statement in this paragraph 
has its counterpart in the next. 

The Maxwell equations for neutral, dielectric media, 
V • D = 0, V B = 0, 

D = V x H M , B = -V x E M , (12) 

impose the analogy g -> D, B and LT y -> H M , E M : The 
thermodynamic variables are D and B, being even and 
odd, respectively. Eqs([l2"|) are their equations of motion, 
while the two non-temporal Maxwell equations are con- 
straints. [Eqs(|l^) must have this form to ensure that the 
two constraints are satisfied at all time. And we already 
know E = E M and H = H M in equilibrium.] The fields 
H M and E M appear only where fluxes do, they therefore 
split into reactive and dissipative parts, 

H M = H + LL°, E M = E + E D . (13) 

The reactive ones are again thermodynamic derivatives, 
while the dissipative fields H D and ~E D are proportional 
to V x E and V x H, respectively, for the same reason as 
above: If these quantities are nonvanishing, the entropy 
is not maximal with respect to D and B, so dissipative 
fields are generated to push D and B toward equilibrium. 
For an isotropic system, we have 

H D = -(«/mo)V x E, (14) 
E D = (/?/e )V x H + 7 VT, (15) 

where 7VT is a cross term, similar to that producing the 
Peltier effect. Again, since (a//Jo)VxE and (/3/eo)VxH 
are of opposite parity under time reversal as D and B, 
respectively, they account for the irreversibility of the 
macroscopic Maxwell equations. The transport coeffi- 
cients are essentially the relaxation time of magnetization 
and polarization, respectively, 

a = tm(1 - Mo/aOj P = tp(1- e /e). (16) 
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[This can be shown with a simple relaxation Ansatz for 
the magnetization and polarization, which to linear or- 
der in lot m and lotp yields the two dissipative terms of 
Eqs(||£|.] 

This concludes the brief presentation of the hydrody- 
namic Maxwell equations, the first side of the complete 
theory^ Before we consider their ramifications, in sec- 



tion 



VII, and discuss the flip side, 



the force on volume 
there 



elements exerted by these fields, in section VIII 
is a prevalent misunderstanding concerning the coarse 
graining procedure that we need to address first. 




VI. COARSE-GRAINING AND 
IRREVERSIBILITY 

A view point one frequently encounters in textbooks 
takes the macroscopic Maxwell equations as the spatially 
averaged version of the microscopic ones. If true, it con- 
tradicts the hydrodynamic Maxwell equations: Compar- 
ing b = -Vxe to B = -Vx E M , it seems compelling 
that (i) E M = (e), and (ii) E M must remain even un- 
der time reversal, as this property cannot be altered by a 
spatial integration. Therefore, E M must not contain any 
odd, dissipative terms such as (/?/eo)V x H. 

Spatially averaging any microscopic equation of mo- 
tion does not usually lead to a coarse-grained, macro- 
scopic dynamics: An initial macrostate contains a large 
number of microstates which generically evolve into fi- 
nal microstates that belong to very different macrostates. 
This lack of uniqueness renders a macroscopic dynamics 
impossible to formulate, since for the prediction of the 
final macrostate one needs the knowledge of the actual 
microstate. [The averaged dynamics is of course unique 
if the microscopic dynamics is strictly linear - but any 
nonlinear term changes this qualitatively. And as empha- 
sized above, the complete microscopic electrodynamics is 
nonlinear.] 

Fortunately, we are only interested in the time evo- 
lution of the field actually measured, not every single 
microscopic one, even if spatially averaged. So we can 
restore uniqueness by taking an ensemble average, the av- 
erage of all microstates contained in a given macrostate. 
Quite frequently, especially if local equilibrium holds, this 
operation alters the time inversion property of the rele- 
vant coarse-grained field, as in the following example due 
to Onsager. 

Consider a macroscopic variable x that vanishes in 
equilibrium. Expand the entropy, S — Sq — 7a; 2 , to quan- 
tify the exponentially diminishing probability of higher 
values of x, and the strongly reduced number of mi- 
crostates compatible with them. Assume that x relaxes 
quasi-stationarily, ie (partial) equilibrium for a given 
value of x is (compared to its slow relaxation time) es- 
tablished instantaneously. Then, given the initial value 
^o(£o)j the vast majority of microstates in this ensemble 
will at to arrive in xq as their excursion peak while un- 
dergoing fluctuations from equilibrium, x = 0, see Fig 2. 



to 

FIG. 2. Furthest point of excursion 

A few will overshoot a moment later to attain higher 
values x > Xq that are far less probable, still others will 
have come from higher values. But the majority will turn 
around at x to move back towards x = 0. 

The ensemble average will follow suit, to arrive at X\ 
just slightly smaller than Xq. Here, the ensemble is con- 
siderably modified by re-establishing equilibrium, as fur- 
ther microstates, with xi as their excursion peak, join in. 
As the new members greatly outnumber the old ones, the 
ensemble average now follows these to move still closer to 
equilibrium. It is this ensemble-averaging at every step, 
over an ever increasing group of microstates represent- 
ing ever more probable macrostates, that increases the 
entropy, forces the macrostate towards equilibrium and 
renders the macroscopic dynamics irreversible. 

In the macroscopic formulation, x = —x/t, this 
physics is accounted for by the dissipative term on the 
right side ~ x. Its form is valid irrespective of the un- 
derlying microscopic dynamics, which certainly does not 
contain any such term. In a sense, x only measures the 
"entropic distance" to equilibrium. 

The same is true for the heat diffusion equation (in 
a stationary medium), where the only current ~ VT is 
again entropic and does not reflect an underlying micro- 
scopic dynamics: The thermodynamic force VT vanishes 
in equilibrium, same as x, and is again a measure of the 
distance from it. So both types of macroscopic dynam- 
ics are in fact quite similar in their construction, with 
the only difference that the variable carrying out diffu- 
sion macroscopically is a conserved quantity. In the cases 
7r *j~r" 7I y ) H + H 15 , E + E D , the respective first term does 
reflect the microscopic dynamics, while the second term 
is entropic. They measure the distance from equilibrium 
and nudge the respective field towards it. 

Summarizing, we conclude that since E M and H M are 
not simply the spatial average of e and h, there is no good 
reason to rule out the dissipative fields E D and H D . A 
corollary result is the difference between the stationary 
solution, V x E M , V x H M = 0, and the equilibrium 
configuration V x E, V x H = 0. The same must hold 
in the linear response theory, but does not in its usual 
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version, as we shall find out in the next section. 



VII. LINEAR RESPONSE REVISITED 

Let us recapitulate the usual derivation of the macro- 
scopic Maxwell equations, due to Lorentz. Start from the 
microscopic Maxwell equations containing only b and e, 
divide the charge into (g e ), V • p, the current into (j e ), 
dtp, V x m, and define h = b — m, d = e + pto eliminate 
p, m while preserving the structure of the original equa- 
tions. These simple, identical algebraic manipulations 
already yield the structure of the macroscopic Maxwell 
equations, albeit in terms of b, h, d, e. Therefore, these 
equations are still detailed, reversible, and completely 
equivalent to the starting point. A follow-up spatial av- 
eraging is easily executed, as the Maxwell equations are 
linear, and substitute (b), (h), (d), (e) for b, h, d, e. Now 
these equations are of low resolution, yet they have re- 
tained the time reversal symmetry. 

The next step, deceptively simple, is the crucial one. It 
introduces irreversibilities and takes us across the Rubi- 
con into coarse-grained, macroscopic physics — although 
this step is usually considered the input of material prop- 
erties, extrinsic to the Maxwell equations proper. We 
identify the variables, B = (b), D = (d), and (instead 
of going through the averaging over ever increasing en- 
sembles) substitute (h), (e) with H M , E M , and take the 
latter as functions of B, D and their temporal derivatives. 
Retaining only first order derivatives, we have 



E 
H 



M 



M 



D/£+G0/e o )D, 
B/fi + (a/no)B, 



or in Fourier space, 

D = eE M /(l - iwpe/so), 
B = p,H M /(l - iwafl/fjio), 



(17) 



(18) 



where e, /2, a, (3 are positive coefficients chosen to coin- 
cide with the hydrodynamic notation. Note especially 
the lack of time reversal symmetry of Eqs(|l7|). 

Two points here need amplification. First, it would not 
have been correct to take D as a function of E M and its 
derivatives: Generally, starting from a temporally nonlo- 
cal relationship between D and E M , both seem possible, 
and one would need microscopic details to decide which is 
realized in a specific case. Furthermore, since the micro- 
scopic field d contains the polarization p, it is considered 
an auxiliary field, and preferred by many to depend on 
the true field e — to lowest order in to in the form 



D = e(E — /3eE /so), 



(19) 



written such that e = D/E M remains unchanged from 
Eq(|l8|) in the given order. However, this formula is un- 
acceptable on general grounds, irrespective of the mi- 
croscopies. Assume homogeneity and stationarity in 



Egs©, V x E M , V x H M , D, B = 0, and find Eq© 
producing a run-away solution E M ~ exp(t/j3e). (One 
must not change the minus sign in Eq([l9|) to render the 
solution decaying, as this would result in self-amplifying 
electromagnetic waves.) 

The second point concerns the lack of spatial nonlo- 
cality, or why the permeabilities e and \x have not been 
taken as functions of the wave vector q. Usually, the an- 
swer entails a discussion of correlation lengths, which (for 
simplicity) takes place in an infinite medium. And the 
result is that spatial nonlocality becomes important only 
at microscopically small scales, so is usually negligible — 
except perhaps in systems with fast moving charge carri- 
ers, such as dilute plasmas. Let us, however, consult the 
hydrodynamic theory the concept and considerations of 
which naturally include boundaries. It is in fact quite 
unambiguous on this point, and it must be completely 
equivalent to the linear response theory in the parameter 
space where both the frequency and the field strength are 
low, and where the two stripes, A and B of Fig 1, over- 
lap. Assuming constant temperature, we may rewrite 
Eqs(pLp|) as 



where 



E D = (/3/e )D+A 2 V x V x E, 
H D = (a/^ )B+A 2 V x V x H, 

A = c\ ad. 



(20) 
(21) 

(22) 



may 



If either a or (3 is very small, the two terms ~ A 2 
be neglected. Then these two equations - via Eqs(||) 
and (H3) - reduce to the usual linear response expres- 
sions, Eqs(|l7],[l8|). However, there are also cases in which 
A cannot be neglected: Being proportional to the relax- 
ation time of magnetization and polarization, cf Eq([l6]) , 
a and (3 vary greatly, from j3 « 10~ 15 s for transparent 
dielectrics, to a w lCP 5 s for colloidal magnetic liquids 
(Rosensweig 1985, Shliomis 1974); while water (or any 
other matrix liquid with a strong, permanent molecu- 
lar dipole moment) is in the middle range, (3 ~ 10 _9 s. 
So a water-based ferrofluid should have a colossal A ~ 
3xl0 3 cm. 

There is a statement about dielectric media, found in 
every textbook on electromagnctism: If the electric and 
magnetic field are static, they are longitudinal and de- 
coupled from each other. It is true if the usual linear re- 
sponse theory Eqs ( |l7| , |T8| ) , holds: Setting D, B = 0, we 
have VxE, VxH = 0, and the only stationary solution 
are indeed longitudinal, decoupled and in equilibrium. 
But it is not true in the general case, as the more com- 
plete version of the linear response theo ry, Eqs(^0| , ^l| ) , or 
simpler, the hydrodynamic theory Eqs ( Jl^ , |l3| , |l 5| ) , enter- 
tains a stationary solution that is off-equilibrium and of 
the form 



E x = {£ + e z t x + £-e- z l x )/^ 
H y = {£+e z ' x - £_ e -*/ A )v/a/(/3 Mo ) 



(23) 
(24) 
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where £± are constant amplitudes. These electric and 
magnetic field are coupled and transverse, they start off 
from the boundary and relax exponentially into the bulk. 
Since the presence of a boundary makes itself felt over the 
distances of A, a dielectric ferrofluid is indeed a dense and 
strongly interacting system that entertains a macroscop- 
ically large spatial nonlocality. It would be interesting to 
detect these fields, and fortunately, this does not appear 
to be difficult. 

[Take a slab of dielectric ferrofluid, with a width L <C 
A, say 1cm. Expose this liquid to an oscillating electric 
or magnetic field, tangential to the slab, of the frequency 
lu <C c/L, and measure the internal field. Conventionally, 
we expect the result to be a uniform internal field, E or 
H, that oscillates in phase with the external one and has 
the same magnitude, so D and B display a phase lag 
~ s",p". The hydrodynamic consideration includes the 
stationary solution above, and the result depends on the 
type of the containing plates: If they are nonconducting, 
the phase lag has the same magnitude but an opposite 
sign; if they are conducting, the internal electric field 
should be drastically reduced (Liu 1997).] 

VIII. ELECTROMAGNETIC FORCES 

Given the input of the thermodynamic theory, the ba- 
sic form of the Maxwell equations and the relevant con- 
servation laws, the derivation of the complete hydrody- 
namic theory is an exercise in cogent deduction and alge- 
braic manipulation (Khalatnikov 1965, Henjes and Liu 
1993). One of the main results is the expression for 
the total stress tensor, Iljj = Kij + n®: The first is the 
Lorentz- Galilean boosted Maxwell tensor, the second is 
the dissipative, off-equilibrium contribution, of the usual 
form and ~ Vij if the liquid is isotropic and the external 
field weak. [Otherwise, it will contain terms ~ VxE, 
V x H, cf Liu (1994).] Inserting the explicit expression 
for Ily into Eq(||), and combining it with and Eqs([l2|), 
we obtain a transparent form of the momentum conser- 
vation, 

P(m + v iVj)[vi + (E x H/c 2 D x B)i/p] 

= ((p e )E+(j e )xB + f D ).-V^. (25) 

The first line contains both the acceleration p%kand the 
Abraham-force J^(E x H/c 2 — D x B). (The latter is 
a small quantity if the electromagnetic wave length of a 
given frequency is large compared to the experimental di- 
mension, as is usual for hydrodynamic frequencies.) The 
four terms of the next line are the proper generalization 
of the pressure gradient and include the reactive pon- 
deromotive forces. They are valid for a two-component 
medium (such as a solution), with p±, pi denoting the 
respective density, and p = p\ + p 2 . Given linear con- 
stitutive relations, assuming that e — eq, p — po are pro- 
portional to one of the densities but independent of the 



other, and neglecting contributions ~ 0(v/c), they re- 
duce to the usual Kelvin force, 

VP — 5 (e — s Q )VE 2 ~\{p- p )VH 2 , 

where P is the pressure in the absence of fields. The third 
line, finally, contains the dissipative stress tensor tt^, the 
Lorentz force — in terms of E rather than E M , in con- 
trast to Eq(^) — and the novel dissipative ponderomotive 
force, 

f D = B x V x H D +D x V x E D . (26) 

It is odd under time reversal, same as 7r^, and accounts 
for forces that arise because the polarization and magne- 
tization are not quite in equilibrium. 

As we shall soon realize, terms in f D of linear order in 
v/c may be just as important as those of zeroth order. 
They have been neglected above to keep the arguments 
and display simple, but are easily retrieved: The fields in 
the Euler Eqs(0), and hence in the dissipative fields, are 
those of the local rest frame, 

H D = -(a/p )V x E°, E D = (f3/e )W x H°, (27) 
E° = E + v x B, H° = H v x D. (28) 

This is plausible: The question whether equilibrium 
is established in a system cannot depend on the ob- 
server's choice of the inertial frame. [The derivation con- 
sists of generalizing the results of Eqs(^,0,|l5|) to systems 
with finite velocities, where a pragmatic combination of 
Galilean and first order Lorentz transformation is em- 
ployed. A fully covariant theory (Kostadt and Liu 1998) 
was also derived to make sure that the additional terms 
are indeed negligible in usual circumstances (Symalla and 
Liu 1998).] 

Assuming for simplicity that (3 (and hence E D ) is neg- 
ligible, and confining our considerations to incompress- 
ible systems of simple geometry — a sphere or a slab such 
that the internal field are uniform — we may rewrite H D 
as 

H D = (a//io)[B-(B.V>], (29) 

where the second term, of order v/c, is a result of the 
fact that H D is given in terms of the restframe field E°, 
rather than E. 

First, consider a solid magnetic sphere, v = fix r, and 
expose it to a rotating magnetic field of the frequency f2 b ■ 
If there are no further perturbations, all forces except f D 
vanish, and H D = (a/po)[(f2B — f2) x B]. The torque 
on the sphere, after being partially integrated (over a 
surface slightly larger than the sphere), is J dV(r x f D ) = 
/ dV(B x H D ). With O denoting the sphere's moment 
of inertia, Eq(||) therefore reduces to, 

& = {aV/Qeo) B x [(f2 B - ft) x B]. (30) 

After an initial surprise, it is reassuring to see how these 
two terms of successive orders in v/c cancel each other 
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if the magnetic field and the sphere co-rotate - as they 
must to conserve total angular momentum. (A curiously 
similar effecr at linearly polarized magnetic field was ob- 
served by Gazeau et al, 1997) 

Next, consider a slab of ferrofluid sustaining a shear 
flow V„W(, where the subscripts denote the normal and 
tangential direction with respect to the slab. If B = 0, 
only the second term of Eq(|29|) remains, which represents 
a correction to the viscosity 77, 

pv t = ?7eV^w t , rje =V + (^/pq)B 2 cos 2 99, (31) 

where tp denotes the angle between B and the slab nor- 
mal. [In the presence of a magnetic field, the viscosity 
77 itself is a function of B and its orientation. Therefore, 
the information about this contribution of a is not easily 
extracted experimentally.] 

If the field oscillates in time, the first term in Eq(ps|) 
becomes important, pv t — r) e V\vt — B n B t V ' n a/ fio- Be- 
cause a (as a function of temperature and density) is 
constant within the ferrofluid, the right side is nonvan- 
ishing only at the surface, where a jumps discontinuously 
to zero. So this force is best accounted for with boundary 
conditions. Solving pv t = rj e V„Vt, the boundary condi- 
tions are 

»7eV„v t = (a/fj, )B n B t) v t =0, (32) 

for a free and sticking surface, respectively. [The free con- 
dition is obtained by setting the off-diagonal stress tensor 
II t „ to zero, while heeding the fact that , rather than 
H t , is continuous across the interface.] 

Expose a slab of ferrofluid with its free surface facing 
upward to a static normal and an oscillating tangential 
field. If the frequency is low enough, we may neglect pvt, 
and the gradient given by Eq(|32|) is constant throughout 
the width L of the slab. Since vt — at the bottom, the 
averaged velocity is (v t ) = (a/ /j,o)B n B t L/2r] e . 

Now forget the gravity, bend this slab into a ring with 
the free surface facing inward, and expose this construct 
to a B-field rotating in the plane of the ring. Take the 
momentary orientation of the field as x, then a counter- 
clockwise rotation has B along y. In the two sections of 
the ring perpendicular to x, the gradient of the velocity, 
according to Eq([j2"|), is positive. Because the velocity 
must vanish at the outer rim of the two sections, it is 
negative in the right section, and positive in the left. 
They combine to yield a clockwise circular flow, opposite 
to the external field. 

In an open, round vessel, the free surface of the fer- 
rofluid (facing upward) curves up at the wall, as wetting 
fluids do. And the capillary region is very similar to such 
a sheet of circular flow, and should rotate opposite to the 
external field — while the bulk of the fluid below rotates, 
more or less, with this region (Rosensweig et al 1990). 

All these effects follow cogently from the identification 
of V x E° as the magnetic thermodynamic force. They 
depend on only one parameter a that has a clear-cut 



physical significance and a definite value. In addition, 
there is a one-to-one correspondence to the analogous ef- 
fects stemming from the dissipative, ponderomotive elec- 
tric force, D x V x E D . 
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